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We present an approach to obtain new forms of the nuclear electromagnetic current, which is 
based on an integral form of the continuity equation. The procedure can be used to restore current 
conservation in model calculations in which the continuity equation is not verified. Besides, it 
provides, as a particular result, the so-called Siegert's form of the nuclear current, first obtained by 
Friar and Fallieros by extending Siegert's theorem to arbitrary values of the momentum transfer. 
The new currents are explicitly conserved and permit a straightforward analysis of their behavior 
at both low and high momentum transfers. The results are illustrated with a simple nuclear model 
which includes a harmonic oscillator mean potential. 
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I. INTRODUCTION 



t^- | One of the oldest and still unsolved problems in nuclear physics is that of the determination of the electromagnetic 
current operator. The basic difficulty lies on the absence of useful constraints one can impose to this operator. For 
a long time, Siegert theorem jlj] has played an important role in this respect, mainly because its application allows 
to avoid the consideration of two-body currents in the calculations. After the application of the theorem, the effects 
' of such currents can be evaluated in terms of the charge contributions, the corrections of which are a priori small. 
However, Siegert theorem only applies in the long-wavelength limit and then it is irrelevant for the values of the 

£ — , ■ momentum transfers usually observed in electron scattering experiments. 

In principle, one can think that current conservation, and more precisely the continuity equation (CE), could 

ON 1 provide the needed procedure. Nevertheless, the task cannot be accomplished due to the impossibility to fix in a 
unique manner the operator, because those terms given as the rotational of any function are not constrained by the 

^ CE - . 

Using current conservation, the longitudinal current is eliminated in terms of the charge operator, and the electron 
r ^ , scattering cross section reads 

^ ', do _ 47T(7 M 

dh ~ f rec {2J, t + 1) 



v L J2 i^a| 2 + v T J2 (i^aI 2 + i^maI 2 ) 

A A 



(1) 



where tc\, tsx, and tM\ are the Coulomb, electric and magnetic multipoles of the transition. It is common to rewrite 
the matrix elements of the electric transverse operators in terms of the charge matrix elements, by using the equations 



U! X /A+l _ 

- q tcx = y 2\ + i tx ~ V ^aTT <a+ (2) 



tEX = i^Ti tx --\i2>rri tx+ - (3) 

Here we have defined the transition matrix elements 

tx± = (Jf\\iTx,\±i\\Ji) (4) 

of the multipole operators 

Tj LM (q) = Jd 3 rj L (qr)Yf L {r)-3(r), (5) 

where J(r) is the electromagnetic current operator of the nucleus and Yj^(f) labels a vector spherical harmonic. 

Eq. (Q) is the multipole expression of the CE. It relates the charge multipoles to the longitudinal multipoles of 
the current. Eq. (0) is the definition of the electric multipoles. Using eq. (0) it is possible to eliminate one of the 
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multipoles t\± in terms of the Coulomb tcx multipole and substitute it in the expression of tsx ■ In this way one can 
use the CE constraint even in the calculation of the transverse form factor. Besides, this procedure allows to minimize 
the errors coming from an insufficient knowledge of the transverse current (e.g. meson-exchange currents) in electron 
scattering calculations. 

This method of rewriting the matrix elements of the electric transverse operator in terms of the charge matrix 
elements, has been applied traditionally in many electron scattering calculations. Despite this ensures that CE is 
satisfied, the procedure is completely meaningless in models which do not satisfy the CE ||. As a consequence, the 
problem of current conservation in model calculations of electron scattering by nuclei is still under theoretical study 

i- §■ . 

A similar situation appears in the relativistic treatment of the off-shell ^NN vertex in (e,e'p) reactions. Several 
prescriptions for the current operators obtained by extrapolating the on-shell currents have been proposed ||. The 
corresponding off-shell vertex operators violate current conservation, which is enforced by elliminating the 3- or 0- 
components using the continuity equation. The problem which arises is that different results are obtained depending 
on the option chosen . 

An alternative approach is provided by the so-called Siegert current (SC) developed by Friar and Fallieros ||[)| by 
extending the Siegert 's theorem to arbitrary wavelengths. This procedure is based on isolating the components of the 
current which are constrained by the CE and replacing them by an adequate combination of the charge multipoles. 
However, the results are not satisfactory [p~0| ] because they show a pathological behavior for high momentum transfer 
when the nuclear model considered does not verify current conservation. This has been also pointed out recently by 
Caparelly and de Passos || in RPA and TDA calculations. 

Two main conclusions of the work of Friar and Fallieros Q must be pointed out. First, the SC appears to be 
preferable to the traditional forms of the nuclear currents even for model which satisfy the CE, because the size of 
the exchange current operators is considerably reduced. Second, they suggest that well behaved current operators, 
satisfying the constraints of current conservation for high momentum transfer, could be determined, but these new 
forms of the current have not been found up to know. 

In this paper we present a new approach to obtain nuclear conserved currents. Specifically we show that the SC 
appears as a particular case of a more general family of current operators which can be obtained from the CE. Though 
all of them produce identical results in models conserving the current, they permit to formulate different prescriptions 
for restoring current conservation in electron scattering model calculations. Besides, our formalism allows us to obtain 
explicitly the asymptotic behavior of these new currents. 

Our method deals with the CE written in integral form (Sec. II) from which the new currents can be obtained 
directly both in coordinate and momentum spaces. Our approach generalizes the Friar and Fallieros procedure and, 
as a consequence, permits to find, a priori, infinite ways of restoring the CE in models which do not conserve the 
current. In Sec. Ill we show how to do it and we study in detail the low and high momentum transfer behavior of 
the conserved currents found as well as that of their multipoles. To finish, in Sec. IV we show some examples in a 
simple nuclear model. The conclusions are drawn in Sec. V. 

II. INTEGRAL CONTINUITY EQUATION 

A. Coordinate space 

Usually, the CE is written in differential form as 

V-J(r) = -i [H,p(r)] , (6) 

where H is the Hamiltonian and p and J are the nuclear charge and current densities, respectively. This equation 
restricts the value of V- J, once the density p and the Hamiltonian H of the system are known, but, as mentioned above, 
it does not permit to fix completely the current because the rotational terms do not contribute to the divergence. 

The key point in our development consists in writing an integral equation equivalent to the CE. To do that, we 
start with the relation 

^[r 2j (r)] = 2J(r) + (r.V)J(r). (7) 
Taking into account the vector identity 

V x [r x J(r)] = [V • J(r)] r - 2J(r) - r • VJ(r) , 
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one can rewrite Eq. (Q) as 

4- k 2 J(r)l = r {[V ■ J(r)l r - V x [r x J(r)]} . 

dr 

Integrating in the radial direction and using the CE (Q) we obtain 

r 2 J(r) - r 2 J(r ) = / dr r {[V • J(r)] r - V x [r x J(r)]} 

Jr 

= - [ drr{i[JJ,/)(r)]r + Vx[rxJ(r)]} 

J ro 

This equation permits to obtain the integral continuity equation (ICE) 

J(r) = lfj(ro)-4 / drr {i[if,p(r)]r + Vx [rxJ(r)]} 

r z r z J rn 



(8) 



(9) 



(10) 



The ICE is equivalent to the CE and if current conservation is fulfilled, the electromagnetic current is recovered by 
computing the right-hand side of the equation. The ICE is simplified by choosing a value of r for which r\ J(r ) = 0, 
such as, for instance, ro = or oo. In addition, the fact that the above expression of the current is manifestly 
conserved, permits to use it as a prescription to restore the CE. This point will be discussed in the next section. 

It is useful to write the ICE in an equivalent form. One starts from the equation 

<l [a 2 J(ar)] = 2aJ(ar) + a 2 r ■ VJ(ar) = aC(ar) , (11) 



do 

where we have introduced the adimensional parameter a and defined the auxiliary current C. It is easy to check that 

C(r) = -i[H,p(r)]r - V x [r x J(r)] 
By integrating Eq. (|ll]) from a — ao to a = 1 we obtain 

J(r) = «oJ(ao r ) + / daaC(ar) 
and introducing the form of C given by Eq. ( |l2| ) we have 



(12) 
(13) 



J(r) = agj(a r) 



H, I daa 2 p(ar) 

J a 



V x 



r x 



daaJ(ar) 



(14) 



As for Eq. ([n]), this equation simplifies by selecting a value of ao such as a§J(aor) = 0. It is of particular interest to 
choose ao = oo. In this case we can write the current as 



where 



Jc(r) 



H, 



J(r) = J c (r)+J m (r), 
da a 2 p(ar) r = i 



dX 
A 1 



p(r/X) 
1 dX. 



J m (r) = V x (r x jf a J{ar)\ = V x X J ^ J ( r / A ) 



(15) 



(16) 



(17) 



These currents coincide with those obtained by Friar and Fallieros (see equations (8a)- (8d) and (4c) in Ref. ||). Then, 
Eq. (|lj) can be considered as an extension of the Siegert's current formulated by these authors. Other values of the 
parameter ao fulfilling the relation agJ(aor) = (such as e.g. ao — 0) provide new currents potentially useful to 
restore the CE in model calculations. On the other hand, the procedure we have followed is considerably simpler. 

To finish this discussion, we remember that if CE is satisfied, the different currents one can obtain from Eq. ([hj 
coincide. Then, for the two particular cases ao = and oo, we have 



J(r) = / daaC(ar) — / daaC(ar), 

Jo J oo 

and, as a result, 

I da aC(ar) = 
Jo 

This is a global consequence of the CE and is verified by any conserved current. 



(18) 



(19) 
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B. Momentum space 



In order to study the behavior of the currents as a function of the momentum transfer q, it is useful to see how the 
equations obtained in the previous subsection read in momentum space. The corresponding version of the ICE can 
be obtained by calculating the Fourier transform in Eq. (fl0|). However we derive it again to illustrate the differences 
between coordinate and momentum spaces in what refers to the procedure followed. 

We start with the equation analogous to Eq. ((7J): 

j-[«J(q)] = J(q)+q-V q J(q). (20) 

If we consider the vector relation: 

V q [q • J(q)] = J(q) + q • V q J(q) + q X V q X J(q) , 

we obtain from Eq. (|2p| ) 

Tj-[«J(q)] - V q [q • J(q)] - q x V q x J(q) , 
and, by inserting here the CE in momentum space 

q-J(q) = [JMq)], ( 21 ) 

we have 

j-[«J(q)] = [H, V q p(q)] - q x V q x J(q) = -C(q) . (22) 

Here we have introduced the auxiliary current C(q) which, as it is easy to check, is the Fourier transform of the 
current C(r) defined in Eq. (Op. The ICE in momentum space is then obtained by integrating Eq. ( p2| ) 

J(q) = f J(<to) - - /\c(q) (23) 

1 

= ^J(qo) + - / dq{[H, V q p(q)] - q x V q x J(q)} . (24) 

9 1 J an 



As in the case of coordinate space, this ICE reduces by choosing a value go f° r which goJ(qo) = 0. 

We can also obtain the ICE in momentum space equivalent to Eq. (|lj). To do this we introduce an adimensional 
parameter A such that 

^ [AJ(Aq)] = -C(Aq) . (25) 

Integrating this equation in the interval [Ao,l], and taking into account the definition of the current C given in 
Eq. (p2|), we obtain for the ICE in momentum space: 



J(q) = A J(A q) - / 1 dAC(Aq) 



(26) 



A J(A q) + 



H, f dA(V qi o)(Aq) -qx / dAA(V q x J)(Aq). (27) 

J An . J An 



It is worth to note that the second and third terms in this equation coincide with the Fourier transforms of the J c 
and J m currents defined in Eqs. (^) and (0). Besides, it can be shown that, for Ao = 0, the corresponding currents 
introduced by Friar and Fallieros || are recovered. 

Finally, the condition analogous to Eq. ( |l9| ) is found by considering the values Ao = and 00 in the ICE (p7f): 

/>oo 

/ dAC(Aq) = 0. (28) 

The integral in this equation allows us to analyze the asymptotic behavior of the currents in momentum space in case 
the CE is not fulfilled. We discuss this point in the next section. 



4 



C. Multipoles of the current 



In this subsection we obtain the corresponding expressions for the multipoles of the current. As it is known [ OJ , 
the electric and magnetic multipoles are linear combinations of the multipole operators defined in eq. (|^). Taking 
into account Eq. ( |l3| ) we have 



TjLM{q) = f d\j L {qr)Yf L {v) • c^J(a r) + f daaC(ar) 
An easy calculation permits to write 



T 



JLM 



(?) 



— Tj LM ( — ) + / da — k Cjlm ( — ") , 
ao \aoJ J aa ol z \aJ 



where we have defined the multipoles of the auxiliary current C as 

CjLM(q) = J d\ 3L (qr)Yf L {v) ■ C(r) 



Finally, writing the integral (30) in terms of the variable A = 1/a we obtain 



TjLAi(q) = A T JLM (Ao(z) + 



dAC. 



JLM 



(29) 



(30) 



(31) 



(32) 



where Ao = 1/ao- Hence, the multipoles of the current J are given by integrating the multipoles Cjlm(i) of the 
auxiliary current C(r). This equation is the one corresponding to Eq. (Eq) verified by the current J(q), and which is, 
therefore, satisfied by each one of its multipoles separately. 
Introducing in Eq. (B2h the variable q' = \q we have finally 



TjLAi(q) = —T JLM (qo) + - ( dq' C,jLM(q') , 

Q Q J q 



(33) 



where we have defined qo = \aq. 

By choosing the values qo = and oo, for which the first term in Eq. ( |3"5| ) vanishes, we obtain a global condition 
for the integral of the multipoles Cjlm 



POO 

/ dqCj LM (q) =0. 
JO 



(34) 



As for the current, this condition is fulfilled when the CE is satisfied. 



III. PRESCRIPTIONS FOR RESTORING THE CONTINUITY EQUATION 

The equations derived in the previous section assume that current conservation is verified. Often, the nuclear 
current is not conserved in electron scattering calculations, where a known density operator p(r), obtained as the sum 
of one-body single- nucleon densities, and a current operator J NC , such as 

V-j( NC >(r) ? -i [H,p(r)} , (35) 

are considered. 

Here we use the ICE to restore the CE. Starting with the non-conserved current jK NC ) and the charge density, we 
calculate the auxiliary current C^ NC ^ as given by Eq. ( p2[ ) and define new currents according to Eq. (^3|). In particular 
we consider the values qo = and oo for which goJ(<lo) = and we obtain 

J<°>(q) = -i /\'C( NC >(q') (36) 

q Jo 

J<°°>(q) = - / d« Z 'C( NC )(q'), (37) 

9 J q 
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where q' = (q'/qjq.. As discussed in subsec. [IB, the current j(°)(q) is the one derived by Friar and Fallieros and 
afterwards used by Friar and Haxton [0 to calculate the electron scattering form factors. 

The second current J(°°)(q) constitutes a new possible prescription in order to restore the CE. Both currents 
coincide only if the model current is conserved; if CE is not verified, both prescriptions are in general different, and 
one has 

j(°°> (q) - J (0) (q) = - / dq' (q') + . (38) 

q Jo 

The asymptotic behavior of the new currents can be determined from their definitions in Eqs. ©-©. It can be 
summarized in the following properties: 

q^O: J(°)(q) ~ -C( NC )(0), J (oo) (q) = 0(1/?) (39) 

q^oo: J(°)(q) = 0(l/ 9 ), j(°°)(q)~C< NC >(q) (40) 



The first condition is equivalent to the Siegert theorem. In fact, from the definition of C(q) in Eq. (22) we have 

-C< NC )(0) = [ff,(V qP )(0)] = i[H,d], (41) 
where d is the electric dipole momentum of the system 



d = J d 3 rrp(r) . (42) 



Therefore the current J*- -* verifies the Siegert theorem because it equals the time derivative of the electric dipole 
momentum of the system. 

The second current, j(°°)(q), however, does not satisfy this theorem. In fact, from its definition, we have 



y* oo 

lim <?J (oo) (q) = / dqC( NC )(q) ^ 
<?^° Jo 



(43) 



Then, as established in (|39j), this current diverges in the origin as 0{l/q) if the original current is not conserved. 
Therefore it should be discarded as a physical current for low momentum transfer. 

On the other hand, a similar non-physical behavior is shown by the Friar and Fallieros current for large 
momentum transfer, because 



/■oo 

lim <zJ (0) (q) = - / dgC( NC %) ^ 0. 



(44) 



This current goes as 0(1/ q) for large q. 

Instead, the current works well in the large q limit. From Eqs. (|36])-(|37|) we can write 



jMfe) = \i rV^C^V) ^ C( NC )(q), 



,, 2 ~ vm, - ™, (45) 

where the change of variables fj, = 1/q, (/ = 1/q' has been used. Then the behavior of for q — > oo is the same as 
that of the auxiliary current C ( - NC \ which shows the adequate one because it is built from the physical charge and 
current density operators and its definition does not include integrations, but just derivatives and products. Its good 
asymptotic properties for large momentum transfer made this current to be a priori useful in this regime in electron 
scattering calculations, without the handicap of the pathological behavior shown by J*- ). The existence of such a 
current was suggested in Ref. pj, but, to the best of our knowledge, this hypothesis has not been proved up to now. 

In order to have a better understanding of the properties of these two currents, we can also obtain the asymptotic 
behavior in coordinate space. A straightforward calculation, similar to that developed for momentum space, permits 
to write: 



r 

r — > oo 



0: J(°)(r) = 0(l/r 2 ), j(°°>(r) ~ J< NC )(r) (46) 
J(0)( r ) _ -C< NC '(r) , J(°°)(r) = 0(l/r 2 ) (47) 



As we can see from these equations, both currents have also the opposite roles. The current is well behaved at 
large distances, but it diverges near the origin. On the other hand, for r — > 0, the current J^ 00 - 1 goes as J^ NC \ which 
is supposed to be well behaved at short distances, but it does not reach zero fast enough for r — » oo. 
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Finally, we summarize the asymptotic properties of the multipoles Tj^(q). By defining the multipoles of the two 
currents J^ ^ and as 

^U) = -~/V4S(s'), (48) 
i f 00 

T%L(q) = -l tfCffifeO, (49) 

y J q 

respectively, we find the following asymptotic behavior for these two sets of multipoles: 

9-0: Tfl M { q )~-Cf ] %{ q ), T%l(q) = 0(l/q) (50) 

?-°o: Tfl M {q)=0{l/q), T^ M (q) ~ c£$(g) . (51) 



These properties are similar to the ones verified by the currents in momentum space given by Eqs. (|3S|)-([10[). 

The multipoles T<°) are well behaved for low momentum transfer, where the Siegert theorem applies, but fail for 
high momentum transfer. As a consequence, the transition matrix elements corresponding to the current 3^ do not 
go to zero fast enough as q goes to oo, because the bound nuclear wave functions have an exponential-like behavior at 
large q values. This is the reason of the pathology observed in this regime in the electric multipoles computed with 
the current developed by Friar and Fallieros Jl^,|| . 

On the other hand, the multipoles T^ 00 ) do not work well near the origin (they diverge as 1/q), but are well behaved 
for high momentum transfer. Hence they can safely be used in this regime in model calculations. 



IV. A SIMPLE MODEL CALCULATION 



In order to illustrate the results quoted in the last section, we have used a very simple model to analyze different 
transitions in 16 and 39 K and which was previously considered in Ref. j|. In this model, the nuclear structure is 
described by means of a single-particle Hamiltonian of the form: 



H 



H 2 r 2 



Vis 1 ■ s . 



2m ' ' ' " ' 2m b 4 

The eigenf unctions, R n i, are harmonic oscillator functions and the corresponding eigenvalues are 



Vc 



j(j + 1) - 1(1 + 1) - - 



The nuclear charge density operator is taken to be the usual one, 



k=l 



while in the nuclear current density operator we have included the well-known convection and spin-magnetization one 
body terms, 



1 1 1 + r| 



jc (r) = y^ 

J M W = E ' 



[5(r-r k )V Th + V Tk 5(r-v k )} 



k=l 



flp- 



1 + rf 



1 - tS 



V x [8{v-v k )a k ] , 



as well as the so-called spin-orbit current, 



(52) 
(53) 

(54) 



k=l 



In the previous equations M k labels the mass of the fc-nucleon, S fe = cr k /2 is its spin and 73 = 1 or — 1 according this 
nucleon being proton or neutron, respectively. Finally, /ip (/in) is the proton (neutron) magnetic moment. 



7 



The model built in this way satisfies the CE. Within it we calculate the electric multipoles corresponding to the 
current 

J( p °)(r) = J c (r) + J M (r) + J LS (r) (55) 

for some transitions in 16 O and 39 K and the results are considered as "pseudo-data". The parameters Vb, Vls and b 
used in the calculations are shown in Table 1 . They were fixed in order to reproduce the energies of the single-particle 
states around the Fermi level in the two double closed-shell nuclei ( 16 and 40 Ca) of interest. 

In what follows we discuss the point relative to the restoration of the CE by using the currents 3^ and 
defined in Eqs. ([56]) and (|37|), respectively. We have simulated the usual situation of a model not verifying the CE 
by eliminating the spin-orbit current from the nuclear current operator: 

j( NC >(r) = J c (r) + J M (r) (56) 

We are interested in analyzing the goodness of the approach sketched in the previous section in retrieving the "true" 
results. 

We focus our attention on the two following electric transitions: (0 + — » 2s 1 / 2 lp 1 /2)i~ m 16 *-* ano - 0-^3/2 ~^ ^ s i/2^2+ 
in 39 K. In particular, we study the electric multipoles 

4"i(?) - \/^i <Jf 11 lT ^- M 11 Ji) ~ vWrr< Jf 11 iT &+i (9) 11 Ji) ' (57) 

where the multipole operators Tjj±\ are given by Eq. (§) and where (k) stands for (PD), (0), (oo) and (NC). 

In Fig. 1 we compare the multipoles for (k) = (0) (full curves) and (oo) (dashed curves) with the results of 

the calculation done with the full model, (n) — (PD) (dashed-dotted curves). The asymptotic behavior discussed in 
the previous section is now apparent. The calculations performed with J' ', that is the Siegert's current of Friar and 
Fallieros j|, are right at low q, but differs notably from the "exact values" in the high q region, in agreement with 
the findings of Refs. JlOlfl. On the contrary, the multipoles corresponding to are wrong for small momentum 

transfer, providing the correct behavior for large q. 

In order to have a better idea of the goodness of both calculations, we show in Fig. 2 (upper panels (a,b) ) the 
quantity 

A<*>(«) = \t%?\q)\-\t%l(q)\ (58) 

for (k) = (0) (solid curves), (oo) (dashed curves) and (NC) (dotted curves). 

As we can see, this last calculation provides the better result except at very low q, where we know it violates the 
Siegert theorem. However, for q > 1 fm _1 is closer to the "pseudo-data" than the other two calculations. 

It is worth to note the situation of the Siegert's current (solid curves). Although this current provides the right 
behavior of the electric multipoles at very low momentum transfer, they show a considerable disagreement with the 
"pseudo data" for q above 1 fin -1 , where the new current we have obtained for (n) — (oo) (dashed curves) gives rise 
to more accurate results. 

In order to go deeper in the analysis, we have done new calculations without considering the magnetization current 
(see lower panels (c,d)). As we know, this current is not affected by the CE and by ignoring it we can test the 
importance of this piece of the current in the results. However, we can see that the situation does not change too 
much and similar comments to those made above can be stated for these new calculations. 



V. CONCLUSIONS 

In this work we have developed a method to obtain new forms of the nuclear electromagnetic current. The approach 
is based on the integral form of the continuity equation and produces, as a particular case, the Siegert's current 
developed by Friar and Fallieros. As in this case, the new currents can be used to restore current conservation in 
those model calculations in which the continuity equation is not fulfilled. Besides, our procedure permit to understand 
in an easy way the asymptotic behavior shown by the currents. 

We have illustrated the method by means of a simple nuclear model based on a harmonic oscillator potential which 
includes a spin-orbit term. The results obtained show that, at least in the cases studied, the multipoles calculated by 
using two of the new currents do not produce better results than those found with a model in which current is not 
conserved. This puts some doubts concerning the procedures of "restoring the continuity equation" . 
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This situation is similar to the one found in relativistic calculations of quasielastic electron scattering by nuclei 
where no tractable approach to treating the off-shell dependence rigorously exists. This makes inevitable the ad 
hoc modifications of the currents in order to recover conservation, but it is not possible to decide which one of the 
different prescriptions is the better. 

In any case, our approach permits other possibilities to define new currents which deserve a more careful analysis. 
Work in this direction is in progress. 



ACKNOWLEDGMENTS 



This work has been supported by the DGICYT (Spain), under contract PB95-1204, and by the Junta de Andalucfa 
(Spain) . 



[1] A.J.F. Siegert, Phys. Rev. 52, 787 (1937) 

[2] J.E. Amaro, B. Ameziane and A.M. Lallena, Phys. Rev. C 53, 1430 (1996) 

[3] E.C. Caparelli and E.J.V. de Passos, J. Phys. G: Nucl. Part. Phys. 23 (1997) 91 

[4] Y. Horikawa and Y. Tanaka, Phys. Lett. B 409 (1997) 1 

[5] H.W.L. Naus, Phys. Rev. C 55 (1997) 1580; Nucl. Phys. A 628 (1998) 275 

[6] T. deForest, Nucl. Phys. A392, 232 (1983); J.J. Kelly, Adv. Nucl. Phys. 23, 75 (1996) 

[7] J. A. Caballero, T.W. Donnelly and G.I. Poulis, Nucl. Phys. A555, 709 (1993) 

[8] J.L.Friar and S. Fallieros, Phys. Lett. B 114, 403 (1982) 

[9] J.L.Friar and S. Fallieros, Phys. Rev. C 29, 1645 (1984) 

[10] J.L.Friar and W.C. Haxton, Phys. Rev. C 31, 2027 (1985) 

[11] T. deForest and J.D. Walecka, Adv. Phys. 15, 57 (1966) 



TABLE I. Parameters of the potential used in the toy 
model considered to discuss the violation of the CE (see text). 



Nucleus 


V [MeV] 


b [fm] 


Vis [MeV] 


16 


-53.6 


1.67 


-4.20 


40 Ca 


-55.7 


1.80 


-1.90 



FIG. 1. Absolute value of the multipoles tEj given 
by Eq. ( |57| ) as a function of the momentum transfer 
q, for the two transition considered in this work. Dot- 
dashed curves correspond to the "exact" model. Full and 
dashed curves have been obtained for the currents j( A °) 
for Ao = and oo, respectively. 

FIG. 2. Differences A given by Eq. (^|) for the two 
transitions considered in this work. Full and dashed 
curves have been obtained for the currents j( A °) for 
Ao = and oo, respectively. Dotted curves correspond 
to the model in which the current is not conserved (that 
is by ignoring the spin-orbit current (^4|)). Upper panels 
(a,b) include the magnetization current ( p3| ) while this 
current has been not included in the calculations plotted 
in the lower panels (c,d). 
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